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Abstract. In the framework of differential Galois theory we treat the classi- 
c's cal spectral problem \P" — it(x)\l/ = A\I/ and its finite-gap potentials as exactly 
^ ) , solvable in quadratures by Picard- Vessiot without involving special functions; 

the ideology goes back to the 1919 works by J. Drach. We show that duality be- 
tween spectral and quadrature approaches is realized through the Weierstrass 
i permutation theorem for a logarithmic Abelian integral. From this standpoint 

we inspect known facts and obtain new ones: an important formula for the 
^-function and 0-function extensions of Picard— Vessiot fields. In particular, 
extensions by Jacobi's ^-functions lead to the (quadrature) algebraically inte- 
grable equations for the ^-functions themselves. 
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1. Introduction 

Initially the method of finite-gap integration was developed in works [48, 29, 
30, 39, 44, 22, 23] as a periodic generalization of the celebrated inverse scattering 
transform method (the soliton theory [46]). In the very first papers on this topic 
[48, 30, 41, 44, 22] it has become clear that analogs and generalizations of the 
soliton potentials are to be smooth and real functions u(x) in the spectral problem 
defined by the Schrodinger equation 

(1) * xx -u{x)y = \V, 

if the continuous spectrum of the problem consists of finitely many 'forbidden gaps' 1 
(lacunae, bands, zones, intervals [44, 39]). This explains the widely used termi- 
nology 'finite-gap', abbreviated further as FG. Ensuing development of the the- 
ory went far beyond equation (1) and took the algebro-geometric characterization 
[25, 37, 38, 24]. Riemann surfaces and their theta-functions have become the 
main subject of study [41, 24, 25]. Appearance of these nontrivial objects is dic- 
tated by the very nontrivial dependence of solution to Eq. (1) upon parameter A 
and the search for this dependence is a starting and prime subject of the spectral 
[44, 29, 30] and algebro-geometric (9-function) [37, 24] approaches. These the- 
ories are referred frequently to as the 0-function integration. Presently, one can 
say that the intense study of equation (1) over the last decades led to the fact 
that its FG-theory has been developed almost exhaustively. In this connection, it 
is, perhaps, not without interest to consider one more view on integration of the 
problem (1). 

1.1. Motivation. In the early 1980s some authors revealed the two old pa- 
pers by J. Drach [20, 21] wherein equation (1) was integrated 'directly' and main 
results of the theory were presented in extremely condensed form. Although these 
works had subsequently received some mention in the literature [26, 17, 46] with 
a special emphasis to the FG-theory ([8, pp. 84-85]; written by Matveev), some 
surprising facet is the fact that more detailed exposition of Drach's ideology has 
not been presented in the modern literature hitherto. The need for such exposition 
is apparent when taken into account that the works 2 [20, 21] themselves contain 
no any explanations or proofs. In this connection, it is pertinent to make up for 
this gap and sketch an appropriate theory. 

The original approach by Drach is to integrate (1) as an ODE. Indeed, equation 
(1) is primarily a differential equation in variable x even though we consider it 
in the algebro-geometric [37, 7] or spectral context [30, 44] in which the ^ is 
viewed, primarily, as function of A. Anyway, in so far as the ^(x; A) is a function 
of two variables, complete theory must explain this duality and therefore provide 
conversion between l x- and 'A-formulae'. On the other hand, integration of linear 
ODEs is the subject matter of the old and well developed differential extension of 
the algebraic Galois theory which is variously known as the Picard-Vessiot theory 
and sometimes as the Lie-Kolchin theory. The main references in this topic are 
monographs [50, 31, 9] and classical works [51, 32]. Strange though it may seem, 
the explicit discussion of a direct linkage between this theory and the modern 



If parameter A belongs to such a lacunae (it is a line segment on the real axe A), then the 
^-function growths unboundedly as a function of x. 

2 After these works Drach had not longer returned to integration of linear ODEs. 
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aspects of integrable models associated to Eq. (1) appeared comparatively recently 
[45, 2]. 

Correlation between the Picard-Vessiot theory and O-function methods brings 
up the following question: what is the relation between these two techniques upon 
applying them to the linear spectral problems, say, (1)? The formal answer (com- 
monly accepted (?)) might be the following: the 8-series (see next section for 
definition) is a special function solving (1). Indeed, an integration theory by Picard- 
Vessiot begins with the precise definition of a used function class [50, 32]. By this 
is meant that, in the strict sense, without definition of the integrability domain, any 
scheme based on use of the formally postulated 8-series should indeed be consid- 
ered as integration in terms of special functions. We shall show that this is not the 
case. For example, the O-function integration of Hamiltonian finite-dimensional 
nonlinear dynamical systems q = V{q) is know to be a manifestation of their al- 
gebraically invariant Liouvillian integrability. On the other hand, integrability of 
such systems is very well known to be related [19] to their representability through 
certain linear equations (Lax pairs): 

{L{q)y = \ty, ^ = A(q)^} =S> L=[A,L] q = V(q). 

The 'nonlinear Hamiltonian' integrability is certain to entail the 'linear Picard- 
Vessiot' one, if only because the logarithmic derivative ^/^> is a rational function of 
dynamical variables. It should be noted here that the development of this ideology 
can serve as a basis for independent concept of integration of spectral problems 
at all. However, we do not touch here on such Hamiltonian systems and algebro- 
geometric (FG) integration of partial differential equations (PDEs). We focus only 
on a linear problem as such, so our main intention with this work is to show that 
the scheme of FG-integration of the linear spectral equation (1) should be separated 
into the two parts: 

(1) The invariant property of equation (1) to be integrable, i. e., Lie-Kolchin's 
solvability of corresponding differential Galois group [32, 50]. 

(2) Representation of differential fields and solutions in terms of those func- 
tions which of inevitably appear in the theory. These are the 8-series. In 
some particular cases the series themselves satisfy the algebraically inte- 
grable Hamiltonian ODEs. 

By the invariance, here and in the subsequent discussion, we shall informally 
mean an independence of representations by theta-functions. Notwithstanding the 
fact that representation of solutions requires introducing the highly nontrivial tran- 
scendental 6-objects, the integrability mechanism itself is very simple. It coincides 
in effect with an elementary solvability in closed form 3 and thereby trivializes un- 
derstanding of the major portion of the FG-theory. In other words, our main 
purpose is to bring the Picard-Vessiot aspects— fields, their extensions, differential 
Galois group, quadrature solvability, etc — into the foreground and, subsequently, to 
get representations for them in the FG-terms— spectral curves, variables 7 fc , 8's, 
etc. The latter objects, to the best of our knowledge, have not received mention 

3 J. Kovacic, in his famous work [35] on p. 4, notes: 'by a "closed-form" solution we mean, 
roughly, one that can be written down by a first-year calculus student'. As we shall see, this 
'definition' is completely compatible with the transcendental 0-function characterization of the 
FG-integration of Eq. (1). Rephrasing, there is a closed form solution (Theorem 4.2) that can be 
verified by a direct substitution into (1) followed by use of the first- year student calculus: algebra 
and differentiation. 
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in the contemporary works on the Picard-Vessiot integration of linear ODEs. See, 
for example, work [49, first sentence in §3(b)], monographs [50, 45, 9], and quite 
voluminous references therein. Partially, some fragments of the theory, in a context 
of the elliptic Lame potentials u = Ap(x), can be found in book [45] and work 
[10]. We note also that FG-potentials are the Abelian functions [41] and Abelian 
extensions of differential fields were already briefly considered by Kolchin himself 
[33]. There is no escape from the mentioning nice applications of the Picard-Vessiot 
theory of Eq. (1) to the supersymmetric quantum mechanics. They appeared com- 
paratively recently and this theme is the subject matter of recent works [1, 2]. 

1.2. Outline of the work. Section 2 contains the background material: a 
sketch of the classical Burchnall-Chaundy theory of commuting operators and its 
modern formulation in the language of theta-functions. 

In Sect. 3 we briefly recall the needed facts from the Picard-Vessiot theory, 
introduce the base differential fields (Novikov's fields), and motivate their hyperel- 
liptic extensions. Then we describe a structure of the differential Galois group for 
FG-potentials. 

Section 4 is devoted to the quadrature (Drach) characterization of the FG- 
integrability and an explanation as to how the known Weierstrass theorem on an 
Abelian logarithmic integral performs the transition and difference between quad- 
rature and spectral mode of getting formulae for the ^-function. 

In Sect. 5 we first recall that the theta- function formulae can be derived from 
quadratures ones [13] and then show the necessity to represent the previous base 
objects in terms of theta-functions and, in particular, to introduce an important 
object — the 1-dimensional section of the theta-function argument with a free param- 
eter. Owing to some differential properties of theta-functions the theory acquires 
very effective form in those cases when jacobians are reducible to a product of 
elliptic curves. 

In Sect. 6 we completely pass to the theta-function representations and give 
an appropriate formulation to the FG-Picard-Vessiot integrability of Eq. (1). This 
provides a nice analogy with solubilities in the simplest integrability domains like 
field C. 

Section 7 explains how the theta-function reformulation of Picard-Vessiot in- 
tegrability transforms into the differential closedness of the theta-functions them- 
selves. This also gives a new treatment to the spectral parameter and a relationship 
of this treatment to the closedness and linearity of some of defining equations. We 
expound results at greater length for the cases when multi-dimensional O reduces 
to the 1-dimensional Jacobian #'s. By way of illustration we exhibit a simplest 
g = 2 non-elliptic potential. 

Differential properties of ^-functions described in previous section allows us to 
take these as a starting point for definition of the functions themselves. It turns out 
that such a view leads again to a Liouvillian extension but the latter is accompanied 
by introducing a meromorphic elliptic integral and brings up some questions about 
differential structures of the multi-dimensional O. All this is expounded in Sect. 8. 

In Sect. 9 we exhibit some counterexamples fitting no to the canonical FG- 
theta-theory but being certainly integrable a la Picard-Vessiot with the solvable 
Galois groups. One of good examples is the famous and fundamental Hermitian 
equation (containing a parameter) very closely related to the theory of Eq. (f ). 

Section 10 contains some conclusive remarks. 
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2. Background 

2.1. Commuting operators. The standard soliton/FG-solutions of integr- 
able equations are known to be defined through the associated linear PDEs for the 
auxiliary ^-function: 

(2) £({£/}; cy* = A*, %=A({u};d x )V, 

where L and A are the ordinary differential (scalar or matrix) operators with coef- 
ficients {[/} being, in general, some functions of variables (x,t): {U} := {u k (x, t)}. 
The set of function {U} is usually termed as the potential. This is because the 
first of Eqs. (2) contains no dt and thereby may be considered as a spectral prob- 
lem defined by the operator expression L. Moreover, the classical (spectral [30]) 
property of the potential {U} to be finite-gap is determined by the spectrum of this 
eigenvalue problem and does not depend on t. On the other hand, the classical 
Burchnall-Chaundy-Baker (algebro-geometric) formulation [16, 7] uses one more 
spectral problem instead of second equation in (2): 

(3) L({U};d x )V = A*, A({u};d x )V = 

In both the formulations the nontrivial theory appears if equations (2) or (3) are 
compatible; there exists a common solution ^ and the potential is subjected to 
certain conditions. These are the well-known commutativities of operators \L, dt — 
A] = for (2) and [£,1] =0 for (3). 

2.1.1. Why one should pass from (2) to (3)? There is a simple explanation 
as to this question. Let the operator A be determined by an hierarchy of some 
integrable PDEs 

(4) u t = K([u\), 

where, as usual in a formal differential calculus [27] and in the sequel, the symbol [U] 
is used to denote the finite set of derivatives {u, U x , U xx , . . .}. Such hierarchies have 
been well tabulated in the literature [19, 28]. Regard Eqs. (2) from the viewpoint 
of their explicit integration (in some sense of the word). In a straightforward 
statement this problem is impossible to solve because t is a hidden variable in 
the first of Eqs. (2). In fact, this variable may be thought of as an additional 
spectral one and the i-dependence of the potential {U} is unknown/undetermined. 
Complexity of the question is not reduced until Eqs. (2) remain partial differential 
equations. Indeed, these equations do not have a general solution expressible in 
terms of any known functions. This is because Eqs. (4), integrable as they are, are 
not solvable in general. The only way, in order to solve the question, is to consider 
some particular situations when PDEs (4) admit transformations into some ODEs. 
Clearly, such a possibility is related to the self-similar reductions of Eqs. (4) and 
the most simple case is of course the reduction to the stationary variable z = x — ct. 
Assuming now the dependence U = u(x — ct), we get, instead of (2), 

L(u(z);d x )* = A*, ¥ t =A(u(z);8 g )%. 

The time t disappears in the first of these equations and therefore there exists 
a solution in form of separability of variables: = T(t) ■ ip(z). Substituting this 



This is so indeed because any additional parameter in coefficients of L{{U~\;d x ) may be 
formally viewed as a spectral variable defining a spectral operator pencil. 
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ansatz into the last equations, we get immediately a separability parameter fj, and 
an exponential dependence T = exp(/it). Equations (2) thus acquire form (3): 

L({U};d z )^ = \ip, (A({U};d z )+cd z )iP = fit/), 

where U ~ u(z). Compatibility conditions of these equations [L, A + cd z ] = 
generate the stationary Lax-Novikov equations F([U)) = [28] and, incidentally, 
different kind reductions may lead to other kind equations. For example, the general 
form of Painleve equations and their l (L, j4)-pairs' can be obtained by this way [47]. 

2.1.2. Algebraic curve and integrals of Lax-Novikov equations. Let n, m be the 
orders of the operators L, A respectively. Assuming that Eqs. (3) are compatible, 
we conclude, by elimination of the "J from (3), that parameters A and /i are related 

by a polynomial dependence W(n, A) = and commuting operators L, A them- 
selves are also tied by the same dependence: W(A, L) is a zero operator [16, 7]. 
The identity L) = implies that its coefficients, being the differential func- 

tions Ej of {[/}, must be free constants Cj independent of x. These yields a set 
of (compatible) ODEs Ej([U]) = Cj providing some integrals of motion Ej for the 
Lax-Novikov equations mentioned above [19]. The further theory requires that the 
potential {U} be the complex analytic function of the complex argument x. 

2.2. Theta-function formulae. Solution to Eqs. (3) is an n-valued function 
of A (and an m-valued function of fi) and this multi-valuedness is related to the 
algebraic equation W^i, A) = 0. This equation, being viewed as an algebraic curve 
over C, defines a compact Riemann surface 7L of a finite genus. It is well known 
that multi-dimensional 0-functions are the universal tool in order to impart a 
'single- valued form' to the analytic apparatus on Riemann surfaces of multi-valued 
algebraic functions [6]. Baker [7], initiated by work [16], was the first to transfer 
this A-multi-valuedness of the W into a single-valued function of a point V G 1Z 
and to construct the function itself through Riemann's O-functions; the very first 
sentence of the work [7] indicates this. Akhiezer [3] arrived at the same objects when 
considering the problem (f ) from completely different — purely spectral — viewpoint. 
In the 1970s all these discoveries were substantially developed, generalized [29, 
30, 41, 23, 37, 38, 25], and acquired their current 0-function form. The recent 
excellent survey by Matveev [43] is, perhaps, the best work both on the history 
of the question and background material. Most general and modern treatment of 
these results goes back to works by Krichever and reads as follows. 

The spectral variable is thought of as the meromorphic function A = X(V) in 
the sense that the complex number A in problem (2) is replaced with an abstract 
coordinate on 7Z: the point V . The variable x is viewed as a parameter now and 
the function as function on 1Z, is the function ^>(V) of an exponential type 
[3, 7] with essential singularities of some prescribed form [41, 37]. The formula 
realization of this result is known presently as a concept of the Baker-Akhiezer 
(BA) function [8]. The operator L as above and its coefficients (the potential {U}) 
is called the finite-gap or algebro- geometric. Of course, we could equally well say 
the same about \i = ii(P) for operator A in problem (3). 



^An independent treatment of the second eigenvalue of the second operator. 
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The structure of solutions is universal [7, 37]. For example, the scalar problem 

- d n d n ~ 2 

(5) X* := — * + u 2 (x) -^- 2 * + ••• + u n {x) * = A * , 

and therefore equation (1), in the class of FG-potentials, has a solution which is 
given, under some normalization, by a general formula: 

Distinctions between different FG-potentials {uk(x)} are only in the changes of the 
associated algebraic curve W(/j,, A) = and the curve itself (its 1Z with a canonical 
homology base (a, b)) determines all the quantities appearing in (6). Namely: 0(2) 
is the canonical g-dimensional theta-series 

(7) e(z) := e(z\U) = J2 e* i{nN ' N)+2 * i{N ' z) 

JVGZ 9 

of g arguments z = (zi, . . . , z g ), built by an a-periods Il-matrix of normalized 
holomorphic Abelian integrals 11(7') = (ili('P), . . . ,il s (7 7 )) on 1Z; symbols like 
(TV, z) denote Euclidian scalar product (N, z) := NjZj; V is a free point on 1Z; 
II('P) is a normalized Abelian integral of the second kind with only pole of the first 
order at a point 'P 00 at which \{Vcc) = 00; the vector U times 2ni is a vector of 
b-periods of this integral; D is an arbitrary constant g- vector. All this terminology 
is exhaustively expounded in the numerous literature (see, e.g., [6, 8, 18, 24, 25, 
28, 41]) and the algebraic dependence W(fi, A) = is referred frequently to as 
the spectral curve. As for the problem (1), this curve constitutes a hyperelliptic 
equation of the form 

(8) ii 2 = (X-E 1 )---(X-E 2g+ i) 

and all the FG-potentials are given by the famous Its-Matveev formula [41, 29] 

d 2 

(9) u(x) = -2— ln6(.Ti7 + D) + const. 

dx z 

Example 1. Most simple and popular example is a 1-gap potential. It is unique 
and is determined by the Weierstrass elliptic curve 

(10) =4(A-e)(A-e')(A-e"); 

we denote its modulus as n = cj'/ui, where u>, u>' are Weierstrassian half-periods 
normalized by the condition SsLI > [4]. Since this case is a 1-dimensional one, we 
replace V H> u and put 6 11(7') = u. Therefore, normalizing integral II, we have 

r n(u + 1) = n(u) 

H(u) = C(2^u)-2 ?? u \ tt 1 

n(u + n) = n(u) + — u = -—, 

\ ltd ZUJ 

where ^(z) := C(z|w,a;') and r\ := 77(0;, w') are the standard objects accompany- 
ing the theory of Weierstrassian function p(z) := p(z\u>,u>') [4, 55]. Putting for 
simplicity D = 0, formulae (6) and (9) become 

u{x) = _ 2 i n e f ±. In) - 2 1 , A) = 9 fe" U i n ) c {<(2 - u) - 2 " u} ^ 



dx 2 \2uj 



8(=l») 



^Jacobian of an elliptic curve is isomorphic to the curve itself. 
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and A = p(2uu) (this is the formula A = X(V) above). This potential is a precise 
equivalent of the classical Lame form u(x) = 2p(x — uj — lu') and 1-dimensional 
0(z|n)-series here coincides exactly with the Jacobi function 6*3(z|n) defined by the 
standard formula (31). 

All the constructions mentioned above— spectral, algebro-geometric, 0-func- 
tion, and their varieties— are completely equivalent [38], which is why we shall refer 
to these approaches merely as spectral for short. 



3. Integrability of equation (1) by Picard Vessiot 

3.1. The function R. The theory of equation (1) is closely related to the 
fundamental linear differential equation 

(11) R xxx -4(u + X)R x -2u x R = 

determining function R = R(x; A). Integrating this equation and denoting an inte- 
gration constant as fx, we get 

(12) t j? = ~RR xx + Ir 2 x + (u + X)R 2 

(see work [13] for a successive derivation of these formulae with use of Lie's sym- 
metries approach.) Then the FG-solutions to the ^-function for Eq. (1), in general 
position fj, ^ 0, are given by the well-known formula 



^(x; A) = y/R(x; A) 



exp 



±fidx 
R(x; A) 



CM [ R x (x;X)±2 t i , 
(13) =CXP 7 2R(x;X) dX - 

Formulae (11)— (13) and their variations have been repeatedly appeared in the lit- 
erature [23, 25, 5, 8, 28]. Their precise meaning, however, lies in the fact that the 
availability of formula (13) itself does not mean any integrability [13]. This is just 
an ansatz for Eq. (1) and its solution should be written down in terms of indefinite 
integrals; otherwise all the procedure would reduce to re-notations. 



3.2. General formula for the 'I'-function. In the language of commuting 
Burchnall-Chaundy operators [16] an explicit formula for the ^-function (including 

(13) ) results from the sequential elimination of derivatives from the pair of 
differential equations (3) down to the formula 

X 

(14) ^ X = G([U];X,^ * = cxpjG([U];X,fi)dx. 

This simple recipe of getting the 'I'-function, perhaps, has no received mention in 
the modern literature [54]; it is, however, implicitly exploited in monograph [28]. 
Elimination of the last derivative leads to equation of the curve W(n, A) = 0. 
The algebraic dependence VK(A, L) = serves, in some of works, as the basis for a 
formal definition of the algebraic integrability [38]. 
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3.3. Novikov's equations and differential field. Motivated by the desire 
to define a differential field over which the integration is performed, we need to know 
the differential structure of function R(x; A), ft is well known that this function is 
a series in A with coefficients being differential polynomials in u(x). Computational 
formulae for these polynomials have been detailed in the famous work [27, formulae 
(8)]. Finite order differential conditions on the potential appear if equation (11) 
has a solution being a polynomial in A [13]. As in the FG-theory, the equation 
(11) is also known in differential Galois theory as the second symmetric power of 
operator (1) [50, §4.3.4], [52, p. 671]. 

Definition 3.1. The potential u(x) is said to be an FG-potential if equation 
(11) has a solution R(x; A) being a polynomial in A. No restrictions on coefficients 
of the curve (12) (8) have been imposed. 

This definition is not a standard one but it is of course equivalent to the spectral 
[30, Theorem 1], algebraic [27, Ch. 3], or quadrature [21, 13] definitions. The only 
exception is a formal O-function (algebro-geometric) setting because it does not use 
spectra, resolvent R, or quadratures. An important point here is the fact that the 
function domain for the potential is not defined as usual in Galois theory but calcu- 
lated. This calculation is a problem of nonlinear integration, that is integration of 
ODEs for the potential u(x). These are the famous Novikov equations [48]. We shall 
call this base field the Novikov differential field A/ \[u\) in u(x) -representation bear- 
ing in mind that u(x) satisfies a Novikov equation F(u,u x ,u xx , . . . ,u x 29+1 ^) = 0. 
The field Af([u]) consists of rational functions of u(x) over C(A) and its derivatives. 
Parameter A, the g constants c k coming from the integral Gel'fand-Dickey recur- 
rence [27], and branch points Ek of the curve (12) belong to a subfield of constants 
for Af([u}). Here is an example of Novikov's equation under g = 2: 

(15) {u xxxxx - \Quu xxx - 2Qu x u xx + 30u 2 u x ) + c x {u xxx - 6uu x ) + c 2 u x = 0. 

Integration constant /i, in an FG-class, is fixed to be dependent on the pa- 
rameter in equation, that is /i = /i(A), and equation (12) turns into the formula 
(8). At the same time, as soon as R(x; A) becomes a polynomial in A it becomes 
a differential polynomial R([u]; A) € Af([u]). It should be noted that one suffices to 
have only one solution of Eq. (11) belonging to Af([u]) (see Sect. 9.1 further below). 

3.4. Picard Vessiot field, constants, and their hyperelliptic exten- 
sion. As usual, the Picard- Vessiot extension 

Afilu])^) := Af([u])($+, ¥*) , 

i.e., the splitting field, results from attaching to the field Af([u]) integrals , J ,± of 
equation (1) and its derivatives [32, 51, 34, 31, 45]. A simplest kind of extensions 
corresponds to solvable cases of the Galois theory and is known as the extension of 
Liouville [50, p. 33], [45, 31]. 

Definition 3.2. An extension Af of the differential field Af is sakito be Liou- 
villian if there exists a tower of fields Af = Afo C Afi C • • • C Af n = Af such that 
A/fc+i = Afk (ipu ) , where d x ip k or d x \mjj k is an algebraic element over A4- 

In other words, Liouvillian extensions are the natural enlargements of the base 
field performed by a step-by-step adjunction of solution to the simplest integrable 
ODEs: the 1st order linear ODEs 



(16) 



i/> x - Atp + B 
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with coefficients (A, B) being algebraic/rational over previous step field. The stan- 
dard adjunctions of an algebraic element a, integral fpdx, or exponential expfpdx, 
where p <E A4, are obtained by putting here (A, B) = (Q,a x ), (A,B) = (0,p), and 
(A, B) = (p, 0) respectively. 

Solvability of equation by quadratures is directly connected with a structure of 
the group constituting a differential version of the polynomial Galois group. 

Definition 3.3. Differential Galois group Gal(A/'( v I / )) of a linear ODE de- 
fined over M is a set of linear transformations of its solutions 'J's that preserve all 
the algebraic (over M) relations among 'J's and their derivatives \&' n ) (differential 
automorphisms group) . 

The Picard-Vessiot extension must have the same set of constants as M{[u\) 
[32], [34, p. 411], [50]. It is known that for equations of the form (1) the Wronskian 
\E ,/ iVE f 2 — is a constant. Taking (13) as one of possible bases for solutions, we 

get 

(17) VtV- - = ^Rl-2RR XX +A{u + \)B? 

= 2/i(A). 

Clearly, it must be a constant of the Picard-Vessiot field whatever the potential 
u(x) may be. Neglecting for the moment other constants, we obtain that field of 
constants C(A) requires adjunction of the constant p: C(A) H> C(A,^). Character 
of this extension is determined by the A-dependence of the function R. It may be 
polynomial, rational, or essentially transcendental. The case of rational polynomial 
R(X) is possible only if its poles do not depend on x\ for we should otherwise have 
the A-dependent poles of the ^-function, which is impossible by virtue of structure 
of Eq. (1). Hence the function R must be an entire function of A. The Galois 
group does depend in general on parameters of equation but we are interested in 
the following cases: 

• When the Lie-Kolchin integrability structure is the same for generic A ? 

Therefore two kinds of theories do exist, according as the field C(A, /i) does not, 
or does, belong to a finite algebraic or infinite extension of C(A). The latter fields 
are excessively general because there are huge varieties of entire transcendental 
functions (without any classification) and they do not produce any restrictions on 
potential. On the other hand, finite extensions are of fixed algebraic (necessarily hy- 
perelliptic) character and calibrated by the only number, namely, by the A-degree g 
of the polynomial R([u]; A). In this case infinite Gel'fand-Dickey recurrences termi- 
nate and equation (12) leads to a finite set of differential restrictions on u{x) in form 
of differential polynomials Ej{[u}) = Cj. This gives in fact yet another indepen- 
dent motivation (a la Picard-Vessiot) for appearance/availability of a polynomial 
solution to Eq. (11). For brevity, we shall adopt however the previous shorter no- 
tation for the base field Af([u]) without explicit indication of its A, /x-dependence 
Af([u]; A, p) or dependencies on other field constants Af([u]; A, p, Cj, . . .), where dots 
stand for remaining constants which arise upon complete integration of a Novikov 
equation. 

3.5. Finite-gap Galois groups. Below is a characterization of the Galois 
group of equation (1) in the class of FG-potentials. Condition on parameter A of 
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being an arbitrary quantity is a fundamental requirement meant throughout the 
paper. 

Theorem 3.4. The Picard-Vessiot extension jV([u])(^ ) is a Liouvillian ex- 
tension of the transcendence degree equal to 1. Associated group Gal(A/"(M)(vI /± )) 7 
under generic A 7^ Ej, is connected and isomorphic to the group (3 = (° Q -i), 
where a G C. For other values of X it is isomorphic to group & — ±1 

Proof. For generic A's integral (13) does not belong to Af([u]). From (13) it 
follows that this extension is Liouvillian. Take the canonical basis of solutions (13). 
Then the three quantities , , are expressed through the transcendent 
being adjoined ^ + as follows: 

(is) = « = ^* + , K-SgH (^0,. 

Clearly, in the case of such A's that // = these relations cease to be valid since 
(\j/+ 5 \j>~) become linearly dependent on each other. The relations should be mod- 
ified and we choose the following basis ^f ± : 

(19) (*-) 2 = i?, = 1. 

In both of these cases we adjoin one transcendental element (respectively): 

X x 

fr+ =e^pf Rx ^ fi dx or M,+ = Vr[^. 



In the latter case the radical yR = can sometimes be element of Af([u]) as an 
example of the Lame equations shows [55, Ch. 23], [40]. 

Let us check invariance of relations (18)— (19) with respect to linear transfor- 
mation of the basis 

ifr- J ^ ^7*+ + (Se- 
this determines the Galois group & = ("5) completely. In case (18) we get the 
following set of equalities: 

a 7 = 0, PS = 0, a5 + /?7 = l, /3 = 0. 

From this it follows that /3 = 0, 7 = 0, and aS = 1. The number a can not be any 
algebraic one; for we should otherwise have a finite Galois group and the algebraic 
^I> ± -solutions to give a contradiction with a single- valuedness of the general formula 
(6). In case (19) we derive that S 2 = 1, 7 = 0, aS — 1, and no conditions on f3 
(except for degenerate cases of curve). Hence in generic case A ^ Ej the group (£» 
is connected. Its possible forms are thus as follows: 

(S = (o a" 1 ) ° r 0= (o A 
where 5 = ±1 (compare with cases 3, 4 in Proposition 2.2 of [45]). ■ 

To summarize briefly, we conclude that independently of the topological genus 
of the curve (8), 'finite-gap' groups Gal(A/'([ii])(vI' ± )) do not depend on parameter 
A and cease to be diagonal and connected only for isolated values of the parameter. 
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COROLLARY 3.5. Equations (1) of the FG-class are factorizable over J\f([u]): 

«--<. + A)-(fl. + i£ ± !)(4-i^S). 

Remark 1. We used nowhere any specific form of the polynomial R. Theo- 
rem 3.4 is easily restated for arbitrary integrable A-pencils of the 2nd order with 
the only condition that R(x; A) £ Af([u]). Recall that the spectral A-pencil is a gen- 
eralization of the canonical spectral equation of the form (5) to more complex (e. g. 
polynomial) dependencies of the differential expression L on the external parameter 
A, that is L({U}; d x ,X) = 0. An example is the well-known spectral A-pencil of the 
form 

(20) ^ _ ^ x - (a 2 - ^A + uv) * = 0; 

u \ u J 

it arises when integrating the integrable nonlinear Schrodinger equation [28, 46, 8]. 
In Sect. 9 we shall consider other examples of the A-pencils (see also [13, Sect. 5(a)]). 

By virtue of structure of the formula (14) Theorem 3.4 has a direct generaliza- 
tion. 

Theorem 3.6. Let Af([U]) be a Novikov field associated with the pair of Burch- 
nall-Chaundy scalar operators (3). Then both of these equations are integrable in 
Liouvillian extensions of the same transcendence degree. Under generic X, /i the 
differential Galois groups of these equations are connected and isomorphic to the 
diagonal groups C5 = Diag(a, j3, ... ,7) C GL(C). 

Novikov's equations are known to be Hamiltonian systems integrable by Liou- 
ville [ ]. However such a way of their integration is not necessary since determi- 
nation of the potential, i.e., construction of the field Af([U]), is given by formulae 
following completely from the 'linear' Picard-Vessiot theory. It does not require 
Haniiltonians. 

We finish this section with a digression to one remarkable example. It is a 
Matveev 1-positon potential given by the seemingly elementary formula [42] 

(21) u = — 21n ra {sin(aa! + b) — ax — c}. 

Surprisingly, in spite of its complete fitting into the integration scheme above, it 
is not amenable to integration by means of any classical algorithm in the Picard- 
Vessiot theory (Kovacic [35], Singer [51]). Indeed, these algorithms are applicable 
to the finite algebraic extensions of C(x), whereas this u S C(x, e^ ax+b \a,c). 

4. Spectral/quadrature duality. An integration procedure 

4.1. Drach Dubrovin equations. The quadrature Drach approach gives a 
very simple explanation as to why and where the fundamental polynomial 

(22) R([u];X) = {X- ll (x))---{X- lg (x)) 

comes from, what its roots j k are, and why these are precisely the quantities that 
complete the indefinite quadratures 7 . In what concerns the spectral viewpoint, a 
remarkable result by Dubrovin [22, 23] is that the quantities 7 fe arise as zeroes of 
a 0-function since they solve the inversion problem of Jacobi [3, 41, 30]. 



An elementary explanation to appearance of these objects (supplemented with Russian 
translation of works [20, 21]) can be found in [12]. 
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Lemma 4.1 (Drach [21], Dubrovin [22]). Functions J k {x) satisfy the system of 
ODEs 



and potential is determined by the trace formula [29] 

3 2g+l 

(24) U = 2^ 7fe (x)- £ B fc . 

fc=i fc=i 

Proof. One inserts (22) into (12) and takes (8) into account. Collecting the 
result in degrees (A — 7 fc )™, one requires identity under arbitrary A. One gets (23) 
and (24). ■ 

Remark 2. We presented such a way of proof because it is applicable to higher 
order operators and even spectral A-pencils. The reason is that the derivation of 
trace formulae is not evident when generalizing. Formulae of such a kind cease 
actually to be the 'trace formulae' because they have no longer Gel'fand's treat- 
ment of the operator trace analogs. They are also not derivable directly from the 
definition of i?-polynomial like (22); illustrative counterexamples can be found in 
work [11]. 

4.2. Weierstrass theorem and new representation for the ^-function. 

Main content of this section was briefly announced in [13] and we shall present 
here the extensive proofs, derivations, and precise correlation between spectral and 
quadrature approaches. 

Theorem 4.2. Letu(x) be an FG-potential corresponding to the arbitrary curve 
(8). Then solution to equation (1) is given by the quadratures 



71(5) y g ( x ) 

w ± /i dz f w ± /1 dz 

: — A w 




(25) ^(xiX)-- 

where w 2 = (z — Ei) • • • (z — E2 g +%) and functions 7 fc = jAx) are determined 
through inversion of the set of indefinite integrals 

9 ? dz 9 ? dz 

(26) V Z 3- 1 — = 2x + a g , ^ z n — = a n+ i, n = 0, 1, . . . , 9 - 2. 

The FG-potential u(x) is determined by formula (24). 

Proof. Let us substitute (22) into (13) and change the integration variable x 
to z. We then obtain the following rules 

-^dx = dlnJJ(A-7 fe ) — » — -dz, 

k 

wherein q\ = (jr—Ei) ■ ■ ■ (7^.— Ezg+x). Abelian integrals of 3rd kind, as appeared in 
(25), result from differential equations (23). Furthermore, substitution of expression 
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(25) into equation (1) leads, to get an identity, to formula (24); this can serve as yet 
another way of derivation of the trace formula. Symmetrizing right hand sides of 
equations (23), we rewrite them down in a form that admits an application of the 
indefinite integration operations, that is (26). This set determines j k as functions 
of x. ■ 

Corollary 4.3. In u(x) -representation the extension Af([u])('i' ± ) requires the 
one quadrature (13). If formula (25) is used then extensions Af([u])( 1 i< ± ) consist in 
an adjunction of a symmetric sum of the logarithmic Abelian integrals. 

Remark 3 (Definition). Special attention must be given to the fact that 
integrability of the 'linear vp' is also algebraic (hyperelliptic) as is nonlinear inte- 
grability of J\f([u]). More precisely, in what follows the term 'algebraic' will mean 
that ultimate answers contain finitely many indefinite integrals of algebraic func- 
tions and inversions of the formers. 

As for transition from primary A-dependence to the x-one and vice versa, the 
duality between spectral and quadrature representations is nontrivial; it is charac- 
terized by the following statement. 

THEOREM 4.4. Quadrature and spectral approaches are equivalent in the sense 
that the explicit transition between formula (25) and its spectral counterpart (see 
below) is realized through the Weierstrass theorem on permutation of limits and 
parameters in a normalized Abelian integral of third kind. 

Proof. Sub-exponential expression in (25) is a sum of Abelian integrals, each 
with logarithmic singularities on 72- at two points: (z,w) = (A, and branch 
place (z,w) — (oo,oo). According to a Weierstrass theorem, we may exchange 
singularities of an elementary logarithmic integral with its limits [6, 18, 41]. In 
our case, these are (z,w) = (7(2), £>(cc)) and a lower limit (z,w) = (a, +/3), where 
2 = (a — Ei) ■ ■ ■ (a — i?2 S +i). More precisely, switching the places 

_^ ( (A, (j,) 
(a,0) J \(oo,oo) 

we obtain that the difference 

70) A 

w + /i dz f ( w + g(x) w + 1 dz 



A w J [ z — j(x) z — a J w 

a 00 

is to be everywhere finite quantity, that is certain holomorphic integral: 

7 (x) 

Ai(X) +A 2 (X)z + --- + A g (X)z9- 1 } — . 

Sum of g such quantities must be a holomorphic integral depending symmetrically 
on 7's. Expression (25) is thus converted to its dual object 

, A A x 

/„„n T + l \ fw±pJx)dz fw±g (x)dz , , , . . 
27) * ± ~ exp - <^ / ^4 — + • • • + / ^-4 — + holomorhic(A, x) } . 

2 [7 z-TxW w J z ~lg\ x ) w 



SPECTRAL/QUADRATURE DUALITY OF FINITE-GAP POTENTIALS 



15 



This is nothing else but the spectral formula by Its & Matveev [30] deserving to 
be mentioned more often. We reproduce 8 their result as it has been written in [30, 
p. 351]: 

A atj 

(28) w(A)= A-ffiLdA, [du J (X) = 0, j = l,...,n, 

J 2y/p(X) J 

Pn Pj-1 

where M (A) = A" + aiA™ -1 + ...+«„ and 

A 



- / 1 a^aTm + A4(A) 



(30) 1>(x,\) =cxp(ia;w(A) + 2^^fe(A) 

^ fc=i 

Meaning of all the quantities presented in (28)-(30) and transition (25) (30) 
are obvious from the context. To put it differently, the permutation theorem of 
Weierstrass, regarding inverse transition (30) — > (25), is a way of doing normaliza- 
tion of periods of these integrals so that all the parameters in the spectral formulae 
(28)-(30) can be 'dumped' to a common multiplication constant for the By this 
means the non-indefinite integrals (27) or (29) with a parametrical dependence on 
transcendental functions 7 fe (x) turn into the indefinite ones (25) of an algebraic 
function. ■ 

Description of invariant property of equation (1) to be integrable has been 
completed and we conclude the section with the comments about fundamental 
difference between spectral and quadrature modes of getting the formulae. 

4.3. On a Riemann surface. At this point not only do analysis on Riemann 
surfaces does not come into play, but also the surfaces themselves do not appear. 
One has just a designation 



1*:= ^(\-Ei)---(\-E 2g+1 ) 

and the theory consists of elementary substitutions (see footnote on p. 3). On the 
other hand, verification of spectral formulae (28)-(30) is a highly nontrivial task 
since they contain a complete set of transcendental objects of Riemann's theory of 
Abelian integrals and differentiation of an integrand containing 7's. In this respect, 
not using the permutation theorem, the 'spectral integral' in (27), that is 

A 

' w ± g(x) dz 



— j(x) w ' 

00 

would be very akin to an integral representation of any special function, say, com- 
plete elliptic Legendre's integral 

1 

dz 



K(x) = 



(i 



y/(l-Z 2 )(l-X 2 Z 2 ) 



We have not found mention of this important result in the literature. See also formula (3) 
in Akhiezer's work [3], 
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The latter is not expressible by means of any finite Liouvillian extension over C(ic) 
since K{x 2 ) satisfies a 2nd order irreducible 2-F1 (5, \\ l|£ 2 )-hypergeometric equa- 
tion [4]. 

Remark 4 (history). Both the integrabilities are due to Liouville but chrono- 
logically, 'linear integrability' (1830-40s) was preceded by more famous nonlinear 
integrability of Hamiltonian systems (1840-50s). Despite the numerous modern 
literature, the fact that these two kinds of integrability are non-casually related to 
the one name Liouville was first observed by Morales- Ruiz [45, pp. 51-52]. 

5. The 0-functons 

By virtue of the fact that extension A/"([u])(^E' ± ) is transcendental, analytic 
representations for the previous formulae require introducing new functions. These 
are the 0-series (7) involved to the theory by Matveev and Its in their famous work 
[30]. 

Theorem 5.1. The Q -function representations (6), (9) are deducible from 
quadrature (25). The expressions (6) and (25) are proportional to each other. 

Proof of this theorem and consecutive derivation of formula (6) from (25) have 
been detailed in [13, §7]. An important point here is the deductive appearance of 
all the aggregates of formula (6) when it is viewed as an axiomatic one: the integral 
II('P), its periods U, and the Abel map -U('P). 

From this theorem, we may draw the conclusion that when recognizing the in- 
tegrability of linear equations the 0-functions themselves are not necessary. They 
realize a step to be considered as the next one after emergence of an integral sym- 
bol / in (13) and (26). 

5.1. Some differential properties of theta-functions. Theorem 5.1 sug- 
gests a search for representation of the field Af([u]) by means of 0-functions. To do 
this require some differential properties of 0-functions and we show further that 
they are available. All the FG-theory tells us that Abelian and BA-functions satisfy 
certain differential identities. By these identities are meant the fact that Abelian 
functions, as theta-function ratios of linear sections of ^-dimensional jacobians of 
curves, have a lot of differential relations between themselves and many of such 
relations have forms of known integrable PDEs [15, 28]. Adding here exponential 
functions of the BA-type, we involve into analysis (L, A)-pairs for these PDEs. 
Moreover, let FG-potential be expressible through the 0-functions of Jacobi. Then 
not only do Abelian and BA-functions satisfy certain differential identities but 0- 
functions themselves also satisfy some ODEs. 

Denote by 0[f] the standard 0-series of Jacobi with characteristics (e, 6) [4]: 

(31) *ffl(x|r) := f c ^+f) 2 ^M*+f)(*+f) 5 

— OO 

i.e., 6[{] = -<9i, 9[h] = 2 , 0[8] = 03, 0[°] = 04- Let # := 0(O|t) be corresponding 
^-constants and 0i(a;|r) stands for ^-derivative of the series — 0[J](x|r). Period of 
the meromorphic elliptic integral is denoted by i] = £(1|1,t). 
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Theorem 5.2. Jacobian functions 0[s\, 0% with arbitrary integral characteris- 
tics are differentially closed over the field of the (d 2 , if) -constants and satisfy the 



autonomous ODEs 

{ dO[f\ e{ 



(32) 



dx 
dx 



9 ±L 

W] 



ir 2 ^} 2 m 2 ■ w - 4 + ^ (^] 4 + ^R] 4 ) } ■ ^ 



where [|] signifies an integer part of the number 8/2. 

These formulae are consequences of more general differential properties of Ja- 
cobi's functions briefly tabulated in [14]. One can see that the similar properties 
are inherent characteristics of the general 0-functions if the g-dimensional jacobian 
is isomorphic to a product of elliptic curves. Many examples of such reductions can 
be found in [8]. 

Example 2. Define a 0-function with characteristics [p] as follows: 



(33) e[g](z|n) := \ {a ' p) e(z + in« + |^|n 



3 7ri(a,z+iilc[) 



Then in the case g = 2 we have the following identity. 
PROPOSITION 5.3. The reduction formula for genus g 



2 under II12 



(34) 



Ma(z+/3+|ar) 



-)em{w\K) + -^ +e 



Differential properties of these 0, & -functions follow completely from Theorem 5.2. 

This example is a rather illustrative one because curves have often symmetries 
and if a genus-2 curve has an involutory symmetry differed from the hyperelliptic 
one (A, fi) <— > (A, — fx) then one can show that its II-matrix is reducible to form (34). 
Formula (34) is perhaps a simplest case of reduction of the theta-functions to the 
two elliptic tori r and >t. More complex equivalent of (34) is presented in [8]. 

Corollary 5.4. Let jacobian of the curve (8) be splittable into a product of 
the elliptic curves and the collective symbols 0, §, r\ stand for arising Jacobi's 
theta-functions and their constants. Then the field ( Cd(9; r d 2 ,77, . . .) is a differential 
extension ofAf([u]) (dots indicate other constants of the field). 

Example 3. Non-elliptic 2-gap potential (9) with = 0[8o] for the reduction 
case (34): 

(35) u= -2\n xx {6 A (Ux + A\t)6 2 {Vx + B\k) -\6 1 (Ux + A\t)6 1 (Vx + B\h)), 

where {r, x, A, B,V} are arbitrary. Since the differential ^-calculus is completely 
at hand, we get a particular but nontrivial example of solution to Dubrovin's ef- 
fectivization formulae for genus g — 2. Recall that the problem consists [24] in 
determination of the 'winding' vector U and is described by a system of equations 
containing the undetermined fourth derivatives of the ©-function [24]. In the ex- 
ample under consideration the ultimate answer turns out to be quite finite but 
somewhat large to display here. Equation for the one sought-for quantity U is an 
algebraic equation of degree 9 (exercise: derive it). 
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5.2. Linearly exponential divisor and 9-representation of Af([u]). Let 
us use notation of formula (6) and plug into (9) an inessential exponential multiplier: 



Introduction of this term is motivated by the fact that solutions for the VP-function 
(6) are expressed not only through the O-functions but involve an exponential 
factor. We shall call the quantity 9(x U+D) e hx , with h, U, and D being constants 
with respect to d, the linearly exponential divisor or ii^-divisor. Let us form a field 
C(£ x ). The following proposition gives a weaker property than Corollary 5.4, but 
it is valid for arbitrary genera. 

PROPOSITION 5.5. The field C(£ x ) is d- differential and finitely generated. 

Proof. Expression (36) satisfies a Novikov equation which has finite order 
2g + 1. Hence the derivatives -£^<d(xU + D) c hx of order n ^ (2g + 1) + 2 are 
expressed rationally through £ x and its lower derivatives. ■ 

The field Qg — Cg(Q(xU + D)) may be considered as a field generated by one 
linear divisor (h = 0). It is obviously that Af([u]) C Qg. Of course, Qg contains 
now not only Abelian functions but this extension is well defined since 0-series 
and b-periods U are computed once u(x) has been given. For this reason, we can 
redefine equation (1) as one given over Qg and thereby we let 



Novikov's fields. The constant A, for the moment, may be disregarded. 

It should be emphasized that in the FG-integration there appears not merely 
an abstract multi-dimensional 0-series but its specification £ x ; the 1-dimensional 
linear section of jacobians. Moreover, the most general 9-function is an object 
defined up to an exponential multiplier (see (33)) so that we may think of it, and 
therefore of the divisor £ x , as a continual generalization of Q-function with discrete 
characteristics. Thus, Liouvillian solutions A/"([u])( v I /± ) are expressed through the 
9 with a linear dependence of its arguments upon x and the 'linearly exponential' 
multiplier 9 e 11 ^-' 31 first explicitly appeared in Akhiezer's work [3] and subsequently 
was axiomatized by Krichever [37]. 

Remark 5. In what follows we shall exhibit that the continually parametric 
object £ x may be introduced in its own right, in a particular case, through the quad- 
rature integrable ODEs. Interestingly, the different kind sections of theta-function 
arguments can lead to other important equations. One remarkable property of such 
a kind appears even in the g = 1 case. Let us consider the function Q(x\t) as a 
function on a simplest (i.e., straight line) section of the 2-dimensional variety {ja- 
cobian <g> moduli space}. Without loss of generality we may impart to this function 
the form 0\(At + B\t). Then this object generates the general Hitchin class of 
solutions to the sixth Painleve equation in a form exactly as does the finite-gap 
formula (9), that is logarithmic derivative of a ratio of entire functions [14]. 



This 'linear exponent' is a result of the contemporary theory; Baker [7] did not specify an 
exponential 0-structure of solutions but just cited to pp. 275, 289 of his [6]. 



(36) 



u{x) 



(ft 

In 9 (x 17 + D) e h x + const. 
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6. Integration as a linearly exponential 0-extension 

In this section we give a formulation of the 0-function scheme as integra- 
tion a la Picard-Vessiot. Let us consider the above Picard-Vessiot extension [50] 
Af([u]) C A/"([u])( v I' ± ) in the representation u(x). This transcendence is Liouvillian 
and dependence of the \E r on parameter A is also transcendental contrary to the 
'rationality' of the field C(A, fi). Meanwhile, based on Theorems 3.4 and 5.1, we see 
that the field has the following structure: 

^>-c.( x + °: n T )) jw -)- 

It immediately follows that if we pass to the 0-representation (37) then integration 
procedure can be reduced to one operation. Namely, a field, over which an equa- 
tion has been defined, is supplemented with an element of the same form as one 
generating the field itself: 

(38) Af([u]) CAA(6) cC d (e(xU + D)e hx ,0(xU + D + ii(T))e u(v)x ). 

It is significant in this viewpoint that problem of the 'linear integration' drops 
out along with the problem of building the base field Af([u]) being treated as a 
problem of the nonlinear integration. In the 0-representation the potential is de- 
termined only by means of operations in the field A/"(0); formula (36). 

Theorem 6.1. All the embeddings 

Af([u]) cN{G) c AA(M)(* ± ) c C d (&(xU + D), Q(xU + D + ii(V)) e^ p)x ) 

are the Liouvillian extensions. 

PROOF. The fact that embedding Af([u]) C 7V(0) is Liouvillian follows directly 
form formula (9). The statement concerning the second embedding is a consequence 
of (13) because R([u}; A) G Af(Q). Rewriting formula (6) in the form 

0(xU + D + il(-p)) e^ v)x = • Q{xU + D), 

we deduce that property for the last embedding to be Liouvillian results from the 
proportionality of (6) and expression (13) for \E' + (Theorems 4.2 and 5.1). ■ 

Theorem 6.2. For generic A integration of equation (1) in the -representa- 
tion is equivalent to a multiplication of an element S(x) generating the field Af(Q) 
by an adjoined linearly exponential divisor: 

PROOF. Consider H(x) = Q(xU + D)' 1 . It is clear that C a (H) = Af(Q). Then 

(39) * ± (x; X(V)) = C± • E(x) ■Q(xU + D± il(V)) e ±n ^^ 

because all the holomorphic/meromorphic integrals on hyperelliptic curves change 
sign under permutation of sheets fj, t— > — /i: we may write ±il(V), ±TL(V) in (39). I 

This theorem has an important treatment: 

• When passing to the -representation the equation (1) is integrated as if it 
were an equation with constant coefficients. Integration procedure is thus 
trivialized under a proper choice of 'domain of rationality'. The inverse 
transform method for the soliton class is a particular case of this general 
construction. 
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Indeed, the simplest FG-case corresponds to the 0-gap one with 7V(G) = C(A) and 
we need only one linear exponent; the construction (39) acquires the form 

(40) A) = C± -S(x) ■ e ±aWx , E(x) = e '* € Af(B) , 

wherein S(x) has the 'same form' as the adjoint element c a ^ x . Adjoining all the 
exponents associated with an 7V-soliton solution (and their varieties like positons 
(21) or rational solitons), we obtain the general 0-gap case. 

Remark 6. The structure of solution (39) in form of simple multiplication 
of elements generating Af(0) and its extension Af([u]){$ ) is not quite typical 
for equations integrable by attaching the linear exponents [35] or, especially, for 
equations with a solvable Galois group [32]. This property owes its origin to the 
availability of A in equation (Theorem 3.4). 

Transition between u- and 0-representations is transcendental with respect to 
A-dependence and other constants of the field. These constants are the II-matrices 
of curves, [/-periods, and vector D. We may therefore trivialize the scheme above 
if we proceed further and redefine equation (1) over (37) as one given over the 
A-pencil (field) of the £3, (^-divisors: 

£ x (p) := e(xU + D + U(V)) c a{v)x . 

Although field Cd( y £, x (T')) is generated by the infinite number of elements, equation 
(1) itself constitutes an infinite A-pencil of equations. Strictly speaking, both spec- 
tral and quadrature approaches require to look upon Eq. (1) as being a differentially- 
algebraic one: differential in x and algebraic (polynomial) in A. Furthermore, by 
virtue of Proposition 5.5, the arbitrary ^(T-^-divisor generates some solution of 
Novikov's equation. It may be fixed by choice of one element of the pencil Xl^'P): 

(41) Af(e) -^C d (0(xU + D + U(V )) e"^5 x ) . 

Having extended the field (41), that is Cg(£ x ('P )) , to the field Cg(£ x (V)), its 
Galois group becomes trivial since the integral of equation (1) is given now in form 
of a ratio of two field elements. 

7. Integrability and differential closedness 

7.1. Differential closure in terms of ^-functions. Attaching the divisor 
£x(V) as transcendental element with a parameter V tells us that it should be 
introduced to the theory as the base function, along with the available 0's without 
parameter V. Owing to Theorem 6.2 this would arrive us at a closed differential 
apparatus (differential closedness) accompanying spectral equation and potential. 
Presently, the general 0-formula realization of this viewpoint is an open problem, 
which is why we illustrate it by cases when g-dimensional 0-function reduces to a 
combination of Jacobian ones. 

At first glance, from Theorem 5.2, it would seem that supplement of the basis 
(32) with -C-r-divisor of the type d\(x — u) e hx requires also adjunction of all the 
functions 6[, 82,3. a(x — u). That no such complication takes place will be apparent 
from the following statement. 
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Theorem 7.1. For Weierstrassian curve (10) with modulus r = lo'/lo one 
fines an elliptic £ X (V)- divisor A by the formula 

L(U|7 

where h is an extra parameter. Then the six functions A, 9[, and 9k (k = 1, 2, 3, 4) 
satisfy the closed autonomous system of ODEs over the field C(rj, $ 2 , #(u)) : 



A(x; u|r) := 9 1 (x — u|t) expf ^t^-t x + hx) , u ^ Zt + Z, 

\0i(UTJ / 



(42) 



d6 k 0' ln n2 9 n 9 m 8fc-28 10fe-28 



9^1 ^ 2 o2 2 1 , f 



i 9A = tt^I g;(u)-e 2 g 3 g4 + e2(u)e3(u)g 4 (u)-e? , 

Adx fli 0i(u)' 0i-(02(u)-0?-0?(u)-02) 

Motivation for the theorem lies in the fact that Abelian integrals, on the one 
hand, are representable in terms of theta-functions and, on the other hand, are 
differentially closed: the base integrals of 1st, 2nd, and 3rd kind form a differential 
basis. Indeed, derivatives of integrals are functions and any function is expressed 
through the two base ones (p, p') (generators of an elliptic functions field) which 
are in turn integrals of exact meromorphic differentials of 2nd kind. However theta- 
function representation for a 3rd kind integral is still lacking in system (32). 

Proof. Meromorphic functions are at hand since they are formed by 6*-quo- 
tients. Derivatives of meromorphic functions are again meromorphic ones but dif- 
ferentiation of 8k(x) generates 9[(x). The quotient 9[/9i is proportional to the 
Weierstrass £- function which in turn represents a meromorphic (i.e., 2nd kind) el- 
liptic integral [4] ; it is alone as genus g = 1. Without loss of generality we may set 
that the missing integral of 3rd kind has two logarithmic singularities. The place of 
one of them may be fixed at x = and the second place can be taken as a parame- 
ter. Call it u. Residues of a corresponding differential are opposite in sign and they 
can be moved to a common multiplication constant. Such an integral HI is unique 
(even for arbitrary g) up to a holomorphic one(s) since all the other integrals are 
expressed through IH = JSL(x;u). Its derivatives are again meromorphic functions 
and the process closes. It will suffice to add an exponent of HI and integral HI itself 
is given by the known formula 

z 

m / x 1 [W + W dz cr(2.T-2u) 2C(2u)x 

(43) m{x;u):=- = In -— — e M ; , 

2 J z — z w a{2x) 

where z = p(2x),w = p'(2x) and z = p(2u),w = p'(2u). Weierstrassian pa- 
rameters (cj,cj'), presented in (43), are replaced by the one quantity r thanks to 
homogeneity relation uj 2 p(uj x\uj, u>') = p(x\l, r) =: p(x\r). Holomorphic integral is 
absent in system (32) but is present in a basis of Abelian integrals. Missing element 
A can be formally adjoined by setting A (x|r) = x and supplementing system (42) 
with equation = 1. Adding to (43) the holomorphic integral hA a and convert- 
ing the right hand side of (43) to the 0-functions, one obtains that adjunction of 
cxp(HI) is equivalent to adjunction of the object A(x;u\t). Differentiating (43) and 
converting it to the #'s, one arrives, after some simplification, at the last equality 
in Eqs. (42). 
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If u £ TLt + Z then the integral HE turns into a logarithm of meromorphic 
function: ln(p(a;) — e). There is nothing to adjoin. ■ 

COROLLARY 7.2. Let potential be a finite-gap one and g-dimensional jacobian 
split into a product of the elliptic curves. Then differentiation of g-dimensional 
functions ®(xU + D) reduces to a set of the '1- dimensional' equations (42). The 
functions 9k, 0[ and A, taken possibly with different moduli t, form the differentially 
closed basis over which every Novikov's equation of order ^ 2g + 1 and problem (1) 
are integrated. 

In the framework of this corollary integrability of Novikov's equations is a 
manifestation of differential closedness of the first two equations in (42). The third 
equation in (42) 'integrates' equations with a parameter and their consequences 
(see examples in [13]). Here, 'integrates' means that all these equations are nothing 
more than combinations of system (42) and its derivatives. Let us consider some 
examples. 

Example 4. The two gap Lame potential 

V xx = {24p(2x|r) + A}*. 

This is a classical example presented in many places [8, 28]. Corresponding solution 
expressed in terms of our objects reads as follows: 

(i4) ^_iAfc±4M h p(2u|T) _ A» + 12* 93 ( T ) 



dx 0i(s|t) ' 3A 2 - 12 2 52 (t) ' 9 v 36A 2 - 12 3 <7 2 (r) ' 

M 2 = (A 2 - 48g 2 (r)) (A 3 - 36g 2 (r) A + 432 53 (r)) ( & (8)) , 

where <? 2 (t), g 3 (r) are the standard modular T-representations for Weierstrassian 
parameters a, b entering into the curve w 2 = 4z 3 ~az — b [4], Novikov's equation 
(15) is satisfied under (c 1; c 2 ) = (0, — 672<? 2 (t)) and i?-polynomial has the form 

R([u];X) = X 2 -^u\+^u 2 -36g 2 (T), 

where u = 2Ap(2x\r). 

7.2. Non-elliptic example. We consider here a potential being no elliptic 
function but expressible through the 'elliptic' ^-functions. 

Example 5. The ^-function for non-elliptic potential (35). It is a nonlinear 
superposition of the one-gap ^-functions. If we denote for brevity z = Ux + A and 
w = Vx + B we then derive that 

n(A)a: 



A(z+ir;ili 




A(w+i;il 2 


\x) -A(«;lli|r)A(«;;il2|x) 


6i(z+ 


hr 


r) e^w+H 


x) - 6Mt)0i{w\*) 



(45) *(x; a) = v \\ " „ e 



where holomorphic integrals ill, il2 are the certain linear combinations iljt = CkjUj 
of the elliptic holomorphic ones ui, u 2 because the curve (8) corresponding to the 
O-function (34) has the form 

(46) /i 2 = A(A - 1)(A - a)(A - 6)(A - ab) =: P 5 (A) 
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and is realized as covers of two tori defined by moduli r and x: 



p(2u 1 |r) + rf*(r) + ^(r) = 3 ( * jffi _^ 3 4 (r) , 
(47) 

P (2u 2 |x) + ^(x) + *&x) = 3 { }~ a)( ;\~ b iU i(x). 



These are the classical formulae by Jacobi [8] presented in terms of Weierstrass' p 
and the pair of branch points {a, 6} and moduli {r, x} are commonly written down 
for one another [8]. All the information concerning this curve can be found in [8] 
and we omit details of some calculations. We need to compute the integral 11(A). 

Abelian integrals for the curve (46) are expressed through 0, Q'\ and therefore 
derivation of the meromorphic integral 11(A) is a routine calculation. An explanation 
is that the reduction of holomorphic integrals to elliptic ones entails the reduction of 
the meromorphic Abelian integrals to the meromorphic elliptic ones. The latter are 
expressed through Weierstrassian £- function, i.e., 8'/8, and meromorphic elliptic 
functions. First translate formulae for cover (47) into the language of ^-functions: 

UH) gjCOfehj = gjCi) flKuik) = #M g|W 

1 1 n(r)ef( Ul \r) (A-o)(A-6)' 0§(t) 9f(ui\r) tf|(x) 0?(u 2 |x) ' 
This is a complete set of equations determining the holomorphic integrals Ui , u 2 as 
functions of A. The point A = 00 corresponds to the values Ui = \ t and u 2 = g^- 
Drop out for the moment indication of modulus t in the following transformation 
of a meromorphic elliptic integral: 

8 = n!(,l + , 3 _3 fl ^!M' 



^4s 3 - g 2 s - g 3 



12 \^ " 4 "" J "^i(ui) 



On the other hand, the 1st formula in (48) supplemented with use of the standard 
quadratic ^-identities brings this integral into the following expression: 



-= canst [s- X -^f b d X^ f / A " °>< A ~» + + V 

J /fMA) 7 ^(l-o)(l-6)A 3^ J V ^(A) 

that is meromorphic integral on the curve (46) with a pole at A = 00 and a surplus 
one at A = 0. Doing the same for the second torus (U2) with modulus x, we 
obtain one more meromorphic integral with the same infinities. Forming their linear 
combination, we can construct the integral with a single singularity at infinite point. 
After some algebraic simplifications the sought-for result becomes: 



i-\r\', 



(49) 



01(U 2 - 




*) 


61 (u 2 - 


\* 





n ( A ) =a '7r? lHT + k'TT? 1 (- + c • Ui + d • u 2 

- ^ (Ul|T) ■■6-^14+ W ( A+ ^ „ +g -u 1+ r-u 2 , 



0i(ui|r) 0i(u 2 |x) A(A-o)(A-6) 

where constants {a, 6, c, <i, p, 5, r} depend only on parameters of the potential (35), 
i. e., on {r, x, A, i3, V}, and are independent of A. Again, after careful co-ordination 
of all the moduli and normalizing constants the direct check of (1), (35), (45), and 
(49) becomes a good exercise in a differential 0-calculus. We mention in passing 
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that this example cannot be elaborated in the framework of the standard elliptic 
soliton theory [8]. 

7.3. A new treatment of the spectral parameter. In proof of Theo- 
rem 7.1 parameter u was the only 'external' parameter independent of 'internal' 
parameters of the curve (moduli). On the other hand, the only parameter being 
external to the field Af([u]) and equation (1) is A. It plays an isolated role. Apart 
from the fact that it is merely present in equation, it may be treated as an object 
arising from the two independent 'mechanisms': 1) adjunction of a transcendental 
element and, on the other hand, 2) differential closedness of all the Abelian inte- 
grals. As evidenced by the foregoing and formula (25), these are the same things: 

• The logarithmic singularity in a canonical integral of third kind is arbitrary 
and independent of moduli. The property of the theory to be integrable is, 
by construction, independent of it. It may therefore always be thought of 
as a {free) spectral variable. 

The converse is also true. Differentiation of the 3rd kind integrals yields other 
integrals and functions. In other words, informally speaking, one may say that 

• Closed class of ODEs integrable through Q(xU + D) is in fact integrable 
in terms of {0, Q' , £ x } and 'owes' to contain an external (except for 
moduli) parameter V . 

Indeed, there is one fundamental logarithmic integral HI for each algebraic curve 
and it has a single parameter V <:=> A. In the elliptic case this HI has form (43) and for 
arbitrary genera it is expressed through the g-dimensional 0's [6]. The one fold log- 
arithmic ^-derivative of (25), i. e., derivatives of IH^'s; V), yields the rational func- 
tions of (z, iw)(7's) and therefore only meromorphic objects remain since integrals 
themselves disappear. As a rough guide we have here d x exp(M) = Vt x ■ cxp(HI). 
This linear in cxp(HI) equation 10 is treated as a spectral one. Coefficients of this 
equation (more precisely its A-independent pieces) can be thought of as potential. 
The quantity V should always be distinguished as an external one because, other- 
wise, the following chain 

external A -<=> ^-closedness -<=>- ''I'-linearity' 

is destroyed altogether. Moreover, the re-dependence <d(xU + D) is not bound to 
be a linear one and spectral equations must not necessarily be of the form (3). 
Counterexamples in Sect. 9 illustrate these points. 

8. Definition of 9 through Liouvillian extension 

Insomuch as the object A(x;u|r) is in fact a theta-function with a parameter, 
we can use the differential equations described above as the basis for a definition of 
the theta's themselves and, in particular, consider character of their integrability. 

Proposition 8.1. System (32) has the two algebraic (rational) integrals 

d\el - $\el = Axdlel, diet - $\el = A 2 $\d\ 

generalizing the famous Jacobi 9-identities when A± ^= 1 ^ A 2 . 

10 An important remark is in order. We might not say the same as applied to the 'pure 
spectral' object (27) since non-indefinite integral remains. Again, Weierstrass' theorem does the 
job. The 'one-fold d x \ which is equivalent here to the 'c^-differential closedness', explains why 
the spectral equations are always linear. 
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Proof. The straightforward calculation shows that d x A\ = d x Ai =0. ■ 

In a nutshell, the differential genesis of the object A is as follows. Function 9[ 
is determined differentially through 9\ . Therefore two functions with two arbitrary 
constants solve the system (32). One of constants serves a homogeneity 9 i-> CO of 
(32). Another one u is non-algebraic and is related to an autonomy of equations 
(32). Hence the two transcendental functions 9\(x — u|r) and 6*2(2; — u|r) remain. 
These functions are represented, up to a shift and holomorphic integral, by the one 
object A(x; u|r). 

THEOREM 8.2. Differential equations (42) are algebraically integrable. 

Proof. By a direct computation one can show that any solution 9k = 9 of the 
system (42) satisfies the same 5th order ODE 

( 5 °) {^(t) ) +8F ^ ' F = Qn6) lx -2 K , -K:=2 V + ±n 2 (di + i)i). 

Therefore, not taking into account u, there is only one essential parameter in equa- 
tions (42), i.e., parameter n. From this it follows that 

F 

f dz 

F = S(.t; a, b, c) : \ — -^^^^=^^^^= = 2ix + c 

J \/ z(z — a)(z — b) 

and integration is thus completed if, according to definition in Remark 3, we adjoin 
the inversion operation S: 

X X 

(51) = cxpJ{j E(y; a, b, c) dy } dx ■ e Kx2+dx+e , 

where a, b, c, d, e are the integration constants. Of course, the inversion function 3 
here bears no relation to ratios of the (9-series. Integration of equations for functions 
9[ and A is obvious. ■ 

Remark 7. In a separate work we shall show that the algebraic integrability 
above can be supplemented with a Hamiltonian formulation X = ttVH(X) to the 
system (42) and its Lagrangian description. 

It is noteworthy, that variable 9 satisfies an equation of fifth order, not third, as 
it would be expected from the well-known p-equation of Weierstrass [4] . Another 
point that should be mentioned here is the fact that algebraic integrability of equa- 
tions (42) leads not merely to the ^-function itself but to the elliptic £ x -divisor and 
even its non-canonical extension by the quadratic exponent. (Notice that constant 
k depends on modulus but constant d is free.) Further, the two-fold integration of 
the transcendental inversion operation in (51) can be reduced to integration of an 
algebraic function — our base operation. 

COROLLARY 8.3. The 9-function can be defined through a meromorphic elliptic 
integral. 

To prove this it will suffice to make the following substitution in formula (51): 

Z{ X \a,b,c) 



1 



E(y;a,b,c)dy = J 



:dz 



yjz{z-a){z-b) 
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By this we obtain somewhat nonstandard way of introduction of the 9- functions. 
To all appearances, Tikhomandritskii [53] was the first to point out a way of defini- 
tion of the 9 through a meromorphic integral 11 but his note [53] went unnoticed in 
the literature. He poses a question about the natural going from elliptic integrals 
to the theta-functions and presents the mode of transition between these transcen- 
dents by introducing the integral of the 2nd kind elliptic integral. Indeed, rewriting 
formula (51) in the following form 



we observe that such a way of introduction of a ^-function is in effect the result 
of Liouvillian extension of a meromorphic integral, i.e., adjoining an exponent of 
integral of such an integral. By this means we may adopt this point as a differ- 
ential definition of the 9 a la Liouville and, subsequently, construct all the other 
objects of the theory: meromorphic (algebraic) functions are the ^-ratios, Abelian 
integrals of 2nd kind are expressed through just introduced meromorphic integral 
(or, which is the same, the 9'), and the 3rd kind integrals are the logarithmic ratios 
of the 0's with free parameters (the A-objects). Holomorphic integrals are of course 
the independent objects; they are not defined/determined through any other ones. 
An important role of a meromorphic integral was already observed by Clebsch & 
Gordan in preface to their book [18, p. VI] on the base of Jacobi's formula 



where Z(u) is a Jacobi zeta-function notation in the theory of elliptic functions [4[. 

The above differential properties of the 0's splittable to 0's raise the question as 
to whether the general multi-dimensional 9-functions admit the similar 'differential 
kind' definition. In particular, whether exists the purely Liouvillian definition of 
an x-section ®(xU + D) like formula (52) 12 or, if any, the closed set of partial 
DEs defining the complete set of the general 9, 0'(z)-functions as ones of the g 
arguments zl Some relations between 0's and meromorphic Abelian integrals can 
be found in lectures by Weierstrass (though no really this point has been mentioned 
in the modern literature) but the question about closed and differentially Liouvillian 
definition (if it exists) of a 4 ff -set of the ^-dimensional 9-functions and associated 
derivatives 9' remains an important open problem. 



^He does not mention the quadratic extension and differential closedness of the set {0 k , £>(}, 
however. 

12 Roughly speaking, one needs an extensive strengthening of Theorem 6.1; whence it follows 

that 



and, since the ^ is an exponent of the 3rd kind Abelian integral (formula (25)), that is 
= exp]II(7's), the object @(xU + D) itself is computed as a 'Liouvillian extension of Abelian 
integrals': 



(52) 




a 








2ln xx O(xU + D) = ln ra * + (ln^*) 2 + const 




which is a reminiscence of formula (51). 
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9. Non- finite-gap integrable counterexamples 

Definition of integrability domain is not a subject of the 0-function techniques. 
Therefore we may generate integrable equations by any way differed from the classi- 
cal FG-structure defined by formula (6) and Theorem 5.1. For example, equations 
coming no from operators L({U};d x ) by taking the canonical spectral equation 
L({U}; dx) 1 ^ = A "J form in general the operator A-pencils L({U}; d x , A) = 0, say, 
(20). Their solution structure is not known a priori 13 . Moreover, we can even 
construct an equation fitting no in the FG-scheme but having the same formal 
O-function form of solution. 

Example 6. Omitting in notation the elliptic modulus r, elucidate the said above 
by the following modification of the 2-gap Lame equation: 

® xx = |24p(2a;) + 8p(2.T-u) + 16p(u)}*. 

It has a solution of the formal 2-gap Baker-Akhiezer form (44): 

(53) tffou)*^), fc = 4£(u)-2C(2u) 

ax ui(x) 

(exercise: check this solution). By Theorems 7.1 and 8.2 this example is alge- 
braically integrable over Ca(p(2x), p(2x — u)) (and over Cq{0i, A), of course) with 
solvable Galois group but it has little in common with commutative Burchnall- 
Chaundy operators, BA-function, or spectral lacunae. Formula (53) shows that 
this ^-function has no even a pole at point 2x = u where potential does. This pole 
depends on a spectral parameter u lying on an elliptic curve. 

Nevertheless this example should not be considered as 'too artificial' because 
there exist the O-function integrable models having algebraic curves with x-de- 
pendent branch-points. Ernst's equations in general relativity [36] provide a nice 
example along these lines. 

9.1. Hermite's operator pencil. Consider now, in the framework of Picard- 
Vessiot theory, Hermitian equation (11) itself. It is not a Burclmall-Chaundy op- 
erator but the operator A-pencil. By virtue of formulae (11)— (13), we define this 
pencil over field Af([u]) and repeat arguments about A-dependence of R. 

One of solutions to this pencil is not exponential but a purely Abelian mero- 
morphic function (22), i.e., differential polynomial 

Ri = R{[u]\\) =: P 3M([u\). 

Since base of solutions to Eq. (11) is {*+, vp 2 .} and g Af([u]), we may 

put the second solution as a square of the BA-function Ri = 4'+, where vfL). is an 
adjoint transcendent 

*+ : = exp J ^^^ dx ' 
The third solution is i?3 = ^ and we obtain that 

P 2 

(54) #i=P, i? 2 = *+, ^3 = tt- 



13 It is well known, however, that Eq. (20) is related to a matrix canonical eigenvalue problem. 
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Hence rationality domain is the same as in Theorem 3.4, that is Af([u])(^> + ) , and 
this extension is a Liouvillian one of the transcendence degree 1 (we consider only 
the generic case A ^ Ej). We therefore can obtain the following result. 

Theorem 9.1. Under the generic A ^ Ej the Galois group of Hermite's equa- 
tion (11) defined over AfQu]) is connected and isomorphic to the diagonal group 
& = Diag(l, a, a' 1 )- Equations (11) is factorizable over field M {[u]) . 

Proof. As in proof of Theorem 3.4 we perform a linear transformation of the 
basis {i?i, i?2, R3} and check invariance of the base differential relations between 
i?'s. Let us take relations of the zero and 1st order in derivatives: 

i?l=P, i?2-R3=P 2 , (-^2)3; = X p — —R2, (-^3)3 — X p R3 1 

which result from properties (54). Using condition fj, 7^ 0, one easily derives that 
admissible transformations are 

Rx^l-Ri, i? 2 ^ a- R 2 , R 3 ^ 5- R 3 

and R2R3 i-> 1 • R^Rd,- Hence 5 = a -1 . Solutions {Ri, R2, R3} are in general 
not algebraic functions, hence a is a free nonzero complex number and we do not 
need further to analyze the remaining relations of second order in derivatives of 
R's (they will be automatically satisfied). This yields a connectivity of the group 
and the matrix 14 Diag(l, a, a' 1 ). Factorizations of Eq. (11) are deducible by use of 
Liouville's scheme since we know solutions to equation. For example 

d xxx - 4(u + \)d x - 2u x = (d x + Px p — \d x {d x - Px p 2 ^ 

These factorizations are not unique because equation has order 3. I 

9.2. Inversion of non-holomorphic integrals. As a last counterexample 
we consider an equation that leads, on the one hand, to a nonstandard case of 
the inversion problem, and, on the other hand, to a nonlinear x-evolution in a 
theta-function argument. As we shall see, there is no essential difference between 
(quadrature) integrability schemes of this example and those of pure FG-potentials. 
This example was already pointed out as non-standard in [54]. 

Example 7. Let us consider the following spectral problem 

(55) * ra = 4*, 

v 

where v = v(x). As long as we have deal with invariant integration of (55) (Sect. 3), 
the theory has just non-essential modifications and we restrict ourselves to writing 
down all of its attributes in a form of references source for the simplest but nontrivial 
case 5=1. We put 

(56) R(x; A) = vX - <j>{x), cf>(x) := 2a(x - b)(x - c), 



All this can also be seen from the fact that Galois group belongs to SL3(C) and transcen- 
dence degree of the extension is unity. 
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where a, 6, c € C. Novikov's equation is the equation v 3 v xxx — &4>v x + 4(f) x v = and 
its integrals are as follows: 

(57) y 2 = X 3 + 3E 2 ([v])X 2 + E 1 ([v])X + a 2 (b-c) 2 , 

3S 2 = ~-vv xx + - v x -2-, Ei = -<pv xx - -cp x v x + — + 2av. 
2 4 v 2 2 v l 

'Trace formula' follows from a direct analogy of (22), i.e., we set R = (A — 7)1;, 
but inversion problem becomes an inversion of the logarithmic integral rather than 
Jacobi problem (26). Indeed, manipulations with integrals E\, E2 show that 

7 

f 1 dz 1 , x — b „ 

( 58 ) / W n ^ = = TTh 7 to + D > 

J z - E 2 ^/\z 3 - g 2 z - g 3 2a[b — c) x-c 

where g 2 = Y2E% - 4E X and g 3 = 4E X E 2 - %E\ - a 2 (b - c) 2 . 

From (58) it follows that the ^-function description undergoes changes since 
x-evolution on jacobian is essentially non-linear. We pass from parameter E 2 to 
g by the rule p{2g) = E 2 and represent the logarithmic integral (58) in terms of 
^-functions of the holomorphic one r: 7 = p(2r). We arrive at a transcendental 
equation determining function r = r(x): 

9i(r + g) Oi{g) a(b — c) x — c 
As a result we obtain that ultimate solution to the ^-function is far from obvious: 



_ Xp 2 (2r) ± y/(x-b)(x-c) A(r;±u) 

* xx - a 2 (x-bf(x-cy W > W 1 ' j ^(2r) Oi(r) ' 

X = p(2u)-p(2g), h = 2 V u 2 . 

A direct check of this solution is a good exercise in theta-calculus. Moreover, 
nonlinearity of this x-evolution is two-fold. Logarithm (58) contains a fraction- 
linear function but the principal nonlinearity comes from a transcendental nonlin- 
earity of r{x). It never becomes linear even though we replace the general case in 
(56), that is 4>(x) = 2a(x — b)(x — c), with the particular one 4>{x) = const. 

10. Concluding remarks 

The properties of ^-functions outlined above differ in a crucial respect from 
classical special functions since the latter ones are defined by ODEs not integrable 
in quadratures over elementary or algebraic functions. For example Bessel's func- 
tions or the Painleve transcendents. Therefore when generalizing rational theory 
(solitons), not only do algebraic integrability takes place for Novikov's equations 
but it also takes place for linear spectral equations and even ^-functions. In all 
these cases the integration procedure has been closed at a single and common step: 
adjunction of the inversion operation H. The elementary theory does not get by 
without inversion as well: 

rational functions --■» In — » inversion --■» exponent --■ » solitons. 

It should be also emphasized that it makes no difference whether 1st kind in- 
tegrals (Jacobi problem) or 2nd, 3rd kind ones have been inverted (see, e. g., (58)). 
The only thing is needed for the (Liouvillian) algebraic integrability: inversion of 
indefinite integrals of any algebraic functions. Roughly speaking, the inversion 
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procedure appearing in 'theta-methods' has also the Liouvillian characterization 
because, according to important Eq. (16), adjunction of any kind Abelian integrals 
above is allowed. Nontrivial examples on inversion of meromorphic integrals can be 
found in monograph [28]; they are associated with the Camassa-Holm hierarchy 
and have also the theta-function description. In the same place quite extensive 
bibliography is presented. In other words, in regard to invariant integrability, the 
choice of the -series or the inversion operation H is just a question of nomencla- 
ture. Introducing is equivalent to removing 7's from formulae like (23)— (25) and 
conversely. As for analytic representation of solutions, the 0-series is of course the 
fundamental object. 
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